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We explore the subleading-iV c corrections to the large-iV c Balitsky-Kovchegov (BK) evolu- 
tion equation by comparing its solution to that of the all-iV c Jalilian-Marian-lancu-McLerran- 
Weigert-Leonidov-Kovner (JIMWLK) equation. In earlier simulations it was observed that the 
difference between the solutions of JIMWLK and BK is unusually small for a quark dipole scat- 
tering amplitude, of the order of 0.1%, which is two orders of magnitude smaller than the naively 
expected 1/^V 2 s» 11%. In this paper we argue that this smallness is not accidental. We provide 
analytical arguments showing that saturation effects and correlator coincidence limits fixed by 
group theory constraints conspire with the particular structure of the dipole kernel to suppress 
subleading-i¥ c corrections reducing the difference between the solutions of JIMWLK and BK to 
0.1%. We solve the JIMWLK equation with improved numerical accuracy and verify that the 
remaining 1/N C corrections, while small, still manage to slow down the rapidity-dependence of 
JIMWLK evolution compared to that of BK. We demonstrate that a truncation of JIMWLK 
evolution in the form of a minimal Gaussian generalization of the BK equation captures some of 
the remaining 1/N C contributions leading to an even better agreement with JIMWLK evolution. 
As the 1/N C corrections to BK include multi-reggeon exchanges one may conclude that the net 
effect of multi-reggeon exchanges on the dipole amplitude is rather small. 



1 Introduction 

Little is known about the features of small- a; evolution in the Color Glass Condensate (CGC) 
picture [1-28] beyond the Balitsky-Kovchegov (BK) truncation [21-25] of the Balitsky hierarchy 
of evolution equations [23-25]. Besides the theoretical work deriving the Jalilian-Marian-Iancu- 
McLcrran-Weigert-Leonidov-Kovner (JIMWLK) equations that summarize the Balitsky hierar- 
chies in a compact form, only a single numerical study of generic properties of the full evolution 
equations is available, carried out by Rummukainen and Weigert [29]. All other studies employ 
some additional approximation, typically in form of the BK truncation or even more schematic 
approximation. The BK truncation, as the Mueller dipole model [3-5] it is based on, explicitly 
neglects 1/N C corrections to the full ln(l/x) evolution of QCD observables at high energy. Never- 
theless, both JIMWLK evolution and its BK truncation correctly reproduce the AVdcpendence of 
the linear Balitsky-Fadin-Kuraev-Lipatov (BFKL) [30, 31] evolution equation in their respective 
low density limits. This implies that in the linear, low density (BFKL) domain subleading 1/N C 
corrections are manifestly absent from JIMWLK evolution. The influence of 1/N C corrections on 
the non-linear part of the full, untruncated evolution equations is much harder to estimate. 
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The only study of the full leading-ln(l/a;) JIMWLK equation available [29] has established, al- 
beit only summarily, that the 1/N C corrections appear to be much smaller than the 1/N% naively 
expected for the gluon-dominated evolution. Instead of expected l/N^ « 10% corrections, the 
JIMWLK solution for the scattering amplitude of a quark dipole on a target nucleus found in [29] 
differs from the solution of the BK equation for the same quantity by only 0.1%. This has es- 
tablished the BK equation as a reasonable tool to predict the energy dependence of CGC cross 
sections, at least after running coupling and some DGLAP corrections are included [32-35]. 1 How- 
ever, the question remains whether the unexpectedly small difference found in [29] is accidental, 
being perhaps due to either some intrinsic properties of the calculated dipole amplitude or to some 
features of the numerical setup used in [29]. In this paper we argue that the smallness of the 1/N C 
corrections found in [29] is not accidental. In fact it is imposed by an interplay of group theoretical 
properties and saturation effects of the CGC. As a result non-linear small- a: evolution turns out to 
be an example of a system in which the 1/N C corrections are much smaller than naively expected. 

We should emphasize that our discussion remains strictly within the context of JIMWLK evolution 
and within that only explores the contextual neighborhood of the BK truncation. The JIMWLK 
evolution equation is valid for scattering on a large target that provides a strong gluon field, e.g. 
for a nucleus with a large atomic number A. It does not include contributions of diagrams which 
are not enhanced by the strong target field or, for a large nucleus, are subleading in powers of 
A. This excludes, right from the start any discussion of pomcron loop contributions, as they are 
not included in the JIMWLK framework. Indeed for small targets with a weaker gluon field like 
a proton, which has A = 1, pomeron loops are not parametrically suppressed anymore. While 
pomcron loop induced fluctuations have also recently been identified in [36] as a source of possible 
large factorization violations for such small targets with some parametric uncertainty, [37] had 
found that running coupling corrections tend to strongly numerically suppress such fluctuations, so 
that we feel that our exclusion of pomcron loops from the analysis should not lead to a very severe 
restriction for the applicability of our results. 

In addition to the large gluon field in the target required for JIMWLK evolution, the BK evolution 
equation induces a correlator factorization assumption that is valid only in the large- N c limit. [See 
Eq. |6]) below.] Hence the BK factorization ([6]) has two types of corrections: those suppressed by 
the powers of A and those suppressed by powers of N c . In this paper we are interested in the second 
kind of corrections only, in 1/N C corrections, which are resummed to all orders in the JIMWLK 
equation but are excluded in the BK equation. 

Even within the purview of JIMWLK evolution we restrict ourselves to a subset of phenomena: We 
only discuss how 1/N C suppressed contributions affect dipole evolution. Quantities that have no 
good approximation in terms of (multi-) dipole projectiles scattering on dense targets are beyond 
the scope of our discussion. An example not addressed here would be pA scattering at high energies: 
any realistic description of a proton projectile lies far outside the standard dipole large N c approx- 
imation, despite the fact that JIMWLK evolution does cover this example faithfully. Obviously, 
in situations like this, where even a leading order large N c dipole description is unavailable any 
discussion of the size of l/iV c -corrections is moot. 

One should also note that JIMWLK and BK equations were both first derived at leading order 
a s ln(l/x), but have a whole tower of a s -suppressed corrections, of which only the next to leading 
order (NLO) terms are partially available. Running coupling corrections [32-34] have been calcu- 

1 Kuokkancn, Rummukaincn, Weigert, in preparation. 
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lated and partial results for the remaining contributions (new physics channels) are available [34,38]. 
These corrections will change quantitative features to some extent, but should not completely dis- 
tort the qualitative structures found at leading order. Our discussion and simulations will therefore 
focus on the leading order situation and only comment on NLO corrections where possible. 

In Sect. [21 we prepare the ground for our arguments, reminding the reader about the differences 
between the BK equation and the JIMWLK equation for a qq dipole scattering amplitude. The 
removal of subleading 1/N C corrections in the BK equation is operationally achieved by factorizing 
the expectation value of the product of a pair of dipole operators into a product of their expectation 
values. For this reason we refer to the BK equation as a factorized truncation of the JIMWLK 
equation. The difference of the unfactorized and the factorized expectation values measures the 
size of the factorization violations. The factorization violation A is defined in Sect. [2] [see l[8|)]. 
where we also present its main features. At one loop accuracy, a vanishing A would imply a 
complete decoupling of all 1 /N c corrections from dipole correlators and is thus the crucial quantity 
to explore. 2 

In Sect.[3]we will clarify the reason for the smallness of the factorization violations observed in [29], 
for "typical" factorization violations A. A more in depth discussion of this issue than that offered 
in [29] must first note that the correlator A measuring the factorization violation itself has in 
fact contributions that do reach all the way up to their natural size of 1/iV 2 in certain regions of 
configuration space. (Note that A depends on three transverse coordinates: the positions of the 
original quark and anti-quark, and of the emitted gluon. Varying those coordinates gives different 
values of A.) However, as observed in [29], the typical contributions to A in the majority of 
configuration space arc in fact tiny compared to 1/-/V 2 . In Sect. [3] we will systematically map out 
configuration space to identify all regions with factorization violations. We will argue on general 
grounds that the factorization violation A is indeed much smaller than 1/N% in the majority of its 
configuration space, in agreement with the result of numerical simulations presented in [29]. We 
will also demonstrate analytically that the evolution kernel wipes out all contributions from the 
only region where the factorization violation A is of the naively expected order l/iV 2 . We thus will 
complete the proof of the statement that 1/N C corrections to BK evolution, which are consistently 
included into JIMWLK evolution, are indeed much smaller than 1/iV 2 . This constitutes our first 
main result. 

The basis of our mapping out of configuration space in a systematic way is the insight that the 
origin of the factorization violations is to be found in a set of group theoretical identities that apply 
to coincidence points of (s-channel) n-point functions involved in the Balitsky hierarchies (i.e., the 
limits in which any pair of the transverse coordinates overlaps) . The identities are shown in (fT2|) . 
They are, by construction, respected in JIMWLK evolution, but automatically broken at the 1/N% 
level by the correlator factorization assumption underlying the BK truncation. 

In Sect. [4] we note that it is possible to extend the BK equations in a minimal manner that reinstates 
these group theoretical constraints for all eikonal correlators in high energy scattering. The inspira- 
tion comes from calculating all involved Wilson line correlators in the quasi-classical approximation 
known as the McLerran-Venugopalan (MV) model [6-8]. One can sum up all Glauber-Mueller 
(GM) multiple rescatterings [39] to calculate various 2- and 3-point functions (see e.g. [12,40-43]). 
Using the resulting correlation functions one can construct the factorization violation A and study 

2 At NLO, running coupling corrections primarily modify the evolution kernel and thus mainly modify how strongly 
a non- vanishing A affects evolution speed [see also Sect. [5]. Other NLO corrections generically introduce new l/N c 
suppressed contributions but are accompanied by an additional power of a s . 
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its properties. This allows us to revisit our earlier general observation on the structure of A in con- 
figuration space and amend it with explicit expressions for the correlators, albeit within a model. 
However, as was noted in [43] and as we will explain in Sect. 14. 1\ one can also insert the 2- and 
3-point correlators obtained in the GM/MV approximation into the JIMWLK evolution equation 
for the 2-point correlator (the lowest order equation in the corresponding Balitsky hierarchy) . One 
can then suggest treating the resulting equation as an evolution equation in its own right [43], 
though no parametric justification/proof of this statement exists. This equation (see (|22j) below) 
is thus only a guess for the evolution equation beyond the leading-^ BK equation. The result 
will be referred to as a Gaussian truncation (GT) of the Balitsky hierarchies or equivalently the 
JIMWLK equation. This Gaussian truncation had been introduced originally in [43] as an "ex- 
ponential parametrization" for qq dipoles and a certain set of other correlators, with an evolution 
equation derived explicitly for the qq dipole operator. On this level it was also explicitly used in [27] 
to unify a diversity of "McLerran-Venugopalan models." 3 

The relationship of the Gaussian truncation to the BK equation turns out to be unexpectedly subtle: 
On the one hand it extends the BK truncation in the sense that it includes a set of subleading 
1/N C corrections, those "minimally" required to reinstate the coincidence limits violated in the BK 
factorization. Consistently, the Gaussian truncation reduces to BK in the large N c limit. On the 
other hand, Eq. (|22j) . the evolution equation in the Gaussian truncation turns out to be equivalent 
to the BK evolution equation with respect to dynamical content. The only changes occur in the 
way this content is mapped onto the expressions for correlators. 

In Sect. 14.21 we compare the factorization violation given by GT and by JIMWLK, and find them 
similar qualitatively, but still quite different quantitatively. Since we view GT as a truncation 
of JIMWLK evolution and hence the Balitsky hierarchies, we will also clarify where GT breaks 
consistency with JIMWLK: GT remains only an approximation to the full JIMWLK evolution. 

Our analytical arguments are complemented in Sect. [5] by a new numerical study of the JIMWLK 
evolution equation that goes beyond that of [29] with simulations on much larger lattices in the 
transverse space, extending the 48 2 — 512 2 -range covered earlier with simulations on 512 2 — 4096 2 
lattices. We emphasize that the simulations presented here are done for fixed coupling only: at 
present the numerical simulation of the exact JIMWLK kernel with the running coupling corrections 
found in [32- 34] would render the numerical cost prohibitive. To efficiently include them would 
require us to find an alternative representation that allows a factorized form of the JIMWLK 
Hamiltonian akin to that used at leading order. This remains beyond the scope of this paper. 
Nevertheless, the additional numerical effort allows us to reduce extrapolation errors considerably 
(they arise mostly from the infinite volume limit as it turns out) and establish reliably that JIMWLK 
evolution is in fact slightly slower than factorized BK evolution: Subleading I /N c corrections indeed 
slow down evolution, just as was observed earlier with running coupling corrections. Evolution 
speed turns out to be particularly sensitive to factorization violations, which is in keeping with the 
integral expressions of Eq. (|3I|) below. At one loop order, we observe numerically a 3-5% slowdown 
induced by factorization violations where our simulations approach the scaling region. We argue 
that running coupling corrections should suppress the UV part of phase space leading to a strong 
reduction of this difference of evolution speeds between JIMWLK and BK. We cannot estimate the 
influence of other NLO corrections which may well have their own offsetting effects, but they should 
not completely distort the leading order picture. We conclude that, while the net effect remains 

3 This is different in content as well as in spirit from the Gaussian approximation discussed in [44], 
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small, 1/N C corrections pull our predictions towards evolution speeds compatible with experiment, 
not in the opposite direction. This qualitative slowdown effect is our second main conclusion. 

In Sect. [5] we concentrate on the physical origin of 1/N C corrections to the BK evolution equation. 
As noted above, the BK truncation reproduces the ^-dependence of the linear BFKL evolution 
equation, corresponding to a two-reggeon state in the t-channel. However, among other 1/N C 
corrections, the BK truncation neglects contributions of multiple reggeon exchanges [45-48]. Some 
of those omitted higher reggeon exchanges, like the odderon contribution corresponding to a C- 
odd three-reggeon exchange [49-53], have been included in the BK-truncated CGC formalism by a 
minimal modification of the truncated evolution equations [54-56] . Higher-order reggeon exchanges 
usually require substantial modification of Mueller's dipole model to a more generic s-channcl 
picture as one generically is required to include 1/N C suppressed multipole correlators on top of 
simple dipoles: see [57] for an analogue of the Bartels-Jaroszewicz-Kwiecinski-Praszalowicz (B JKP) 
evolution equation [45-48] in the s-channel formalism. Generically not much has been done to 
identify the contributions of higher n-reggeon exchange contributions [58] to nonlinear JIMWLK 
evolution in any systematic way. 

Nevertheless, the fact that the odderon [54-56] and 4-reggeon [57] exchanges are included in the 
s-channel evolution picture allows us to conjecture that all multi-reggeon exchanges are included 
in the JIMWLK evolution equation. The JIMWLK equation also probably includes some multi- 
reggeon vertices containing more legs on the target side of the evolution than on the projectile side. 
If this conjecture is true, one concludes that the difference between the dipole amplitude given by 
BK and by JIMWLK is at least partially due to an aggregate of multiple-reggeon effects. The 
smallness of this difference then, in turn, would indicate the smallness of multiple-reggeon exchange 
effects. 

The link of multi-reggeon exchanges with subleading l/N c corrections gives a natural explanation 
for the slowdown of JIMWLK evolution compared to BK observed in Sect. [5] Generically one would 
argue that nonlinear effects will work to temper any influence of multi-reggeon contributions, which 
would complement the power suppression of 1/N C contributions via the kernel observed earlier in 
our line of argument. If true, this is testable numerically, but it is not a priori clear how to test 
this. Identifying the Gaussian truncation with iterated two reggeon exchange gives a handle on 
this as well: we may filter out the multi-reggeon exchanges by comparing the Gaussian truncation 
with full JIMWLK evolution. It turns out that the Gaussian truncation has a distinctive feature 
that is naturally violated by multi-reggeon exchanges: the Gaussian truncation would predict 
strict Casimir scaling of dipole correlators in different representations. (Casimir scaling is defined 
in (|23p .) In Sect. [5] we illustrate this statement by extending GT to include the simplest multi- 
reggeon contribution in the form of an odderon exchange: we then show that it indeed violates 
the Casimir scaling. Therefore we argue that the size of Casimir scaling violations can quantify 
the net contribution of all multi-reggeon exchanges. We thus can numerically explore the effect of 
multi-reggeon exchanges by measuring the violations of Casimir scaling of the dipole correlators. 
Casimir scaling violation in the numerical solution of JIMWLK that we performed is studied in 
Sect. [6] It turns out that the Casimir scaling violations (which summarize the collective effect of 
all multi-reggeon exchanges included in JIMWLK evolution) are generically small and do not grow 
with energy (see e.g. Fig. I lOj) . This is our third main result. 

We review our results and methods in Sect. 
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2 Dipole evolution in JIMWLK and BK frameworks 



JIMWLK evolution is equivalent to sets of coupled infinite hierarchies of evolution equations, the 
simplest of which is based on the equation for the gg-dipole correlators (S^ y )(Y) for the scattering 
on a target at high energies in which the scattering of the q and q is expressed via light-like Wilson 
lines in the fundamental representation U x or U y respectively (at fixed transverse positions x,y), 



tr^t/t) 

Nr. 



s% ■= \ X T VJ ■ (1) 



This operator is gauge invariant in the sense that the contributions that close the trace at x + = ±00 
are unity to leading order in ln(l/a;). 4 Averaging the operator in Eq. ([1]) over all states in the target 
wave function yields the ^-dependent S-matrix for the scattering of a dipole on that specific target. 
The evolution equation for this average, (S xy }(Y), involves a gluon Wilson line operator U in the 
adjoint representation on its right-hand side. At fixed coupling can be written either as [18,23,24] 

±(tr(U x Ul))(Y) = ^ J d 2 zJC xzy ({[U z ] a \r(t a U x t b Ul))(Y)-C f (tr^ (2) 

or, using |T]) and the Fierz identity 

[U z ] ab 2ti-(t a U x t b Ul) = tr{U x Ut) tr(U»Ul) - ^tr([/ x C/t) (3) 

as 

Jy(S*v)( Y ) = ^Jd 2 z IC xzy (S xz S zy S xy )(Y) . (4) 

The integral kernel in both ([2]) and ([4]) is given by [3,21] 
(x-yf 

xzy - {x-zf {z-yf 1 ' 

Eqs. @ and ^ arc completely equivalent versions of the first equation in the Balitsky hierarchy of 
the quark dipole operator {l}. Eqs. ([2]) and ((4]) obviously do not represent closed equations since the 
evolution of (tT(U x lI^))(Y) depends on an operator with an additional gluon operator U insertion. 

The evolution equation of that new operator, ([U z ] ab tr(t a U x t b U^)) (Y) , in turn will involve yet one 
more insertion of a gluon operator U, iteratively creating an infinite coupled hierarchy of evolution 
equations, the Balitsky hierarchy of the quark dipole operator |T]) [23,24]. JIMWLK evolution 
summarizes the totality of all such hierarchies, based on any (gauge invariant) combination of 
multipole operators but can only be solved numerically [29] at considerable numerical cost. The 
situation can be simplified for the price of introducing an additional approximation that truncates 
the hierarchy. The most widely used truncation is known as the BK approximation. It assumes the 
factorization 

{S xz S zy ){Y) - (S XZ ){Y) (S zy )(Y) , (6) 



4 This is true strictly speaking only if one does not force the {/-factors shown to be one by gauge choice. In this 
case one would need to display the contributions that connect the trace at 2+ = ±00 to a closed Wilson loop. They 
would carry the full contribution, but would remain independent of the paths used to connect x to y. 
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which turns Eq. (fj| into a closed equation in terms of (S xy )(Y) only and thus decouples the rest of 
the Balitsky hierarchy. The BK truncation is valid and is parametrically justified in the large-./V c 
limit for scattering on a large dilute nuclear target. Using j6]) in ((4| we obtain the BK evolution 
equation 



d 2 z K, 



xzy 



(S XZ )(Y) (S zy )(Y) - (S xy )(Y) 



(7) 



Provided that the dipole correlator shapes {S xy )(Y) are not too different (this notion with be 
refined in Sec. [5]) in JIMWLK (without factorization ([6|)) and BK (with factorization, as shown in 
(O) the factorization violation that creates the difference between the two 



A xzy (Y) := (S xz S zy )(Y) - (S XZ )(Y) (S zy )(Y) 



(8a) 



can be simply interpreted as the difference of the correlators on the right-hand side of Eq. ((2]) and 
its BK counterpart ([7|), i.e. 



A xzy (Y) = [(S xz S zy - S xy )(Y)} - [{S XZ )(Y) (S zy )(Y) - {S xy )(Y)} 



(8b) 



(the (S xy )(Y) term is the same under these conditions and cancels trivially) or as a fluctuation 
away from a mean field value 



& xzy (Y) = ((S xz - {S xz )(Y))(S zy - (S zy )(Y)))(Y) 



(8c) 



The first interpretation directly leads us to consider factorization violations as a source for a differ- 
ence in evolution speed in JIMWLK and BK, the second interpretation will lead us to the question 
of what kind of degrees of freedom (which are absent in BK but included in JIMWLK) would be 
associated with these fluctuations. We meet the latter question repeatedly in all remaining sections, 
here we will first look at the individual terms in Eq. (|8bj) to get a generic idea of the structure of 
configuration space and how it affects evolution and then give a first glimpse at how JIMWLK 
evolution via ((4]) might differ from BK evolution (|7|). 

In both cases in an otherwise translationally invariant system with a given parent dipole the in- 
tegrands (correlators and kernel separately - that is why we will leave the latter aside) have a 
twofold mirror symmetry in the z-planc: one with respect to the (x — y)-axis, the other with re- 
spect to an axis perpendicular to x — y, through the midpoint (x + y)/2. The latter only holds if 
(tr(U x Ul))(Y) = (tr([/ x [/t)t)(y) = (tr(U y Ul))(Y), i.e., it is real (and thus symmetric in x <-> y) 
as is the case if we study its contribution to the total DIS cross section at high energy [21,23]. In 
this context it is useful to introduce a z coordinate with respect to (x + y)/2 as the origin 



(x + y)/2 . 



(9) 



There are strong zeroes in the correlators on the right-hand side of the evolution equations as well 
as in A when z — * x or y. They are needed to cancel the kernel singularities at these points and 
have their origin in real virtual cancellations. Generically, these zeroes are not isolated but lie on 
lines that separate the positive from the negative contributions to the integrand of the evolution 
equations. Picking out the positive sign regions in the integrand, the mirror symmetries allow two 
situations: one in which there are two separate such regions adjacent to the q and q respectively, 
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and another where the regions are joined together [generically in situations with dipole correlators 
not too dissimilar from the Golec-Biernat-Wusthoff (GB-W) case [59] (also known as Glauber- 
Mueller multiple rescatterings) which serves as our initial conditions; the initial conditions are in 
fact radially symmetric]. The generic patterns are shown in Fig. [TJ which presents contour plots 
of the right-hand side of the BK equation ([7]) (divided by (S xy )(Y)) obtained by performing a 
numerical solution of that equation. The horizontal and vertical axis on each panel show z' x and 
z' 2 , which are the two components of the two-dimensional vector z'. The coordinates are plotted 
in the units of the initial correlation length R S (Y ) of the system (formally defined as the distance 
at which the dipole correlator falls to 1/2, i.e. via (S\ r \ = R 3 (Y)){Y) = 1/2). R S (Y) can be thought 
of as the inverse of the saturation scale Q S (Y): R S (Y) = 1/Q S (Y). Y is the rapidity of the initial 
conditions for ([7]). The dots in Fig. [JJ denote the positions x and y of the quark and the anti-quark, 
with \x — y\ taken in Fig. [TJ to be equal to R s (Yq). They lie on the contour lines that separate 
positive from negative regions. The left panel of Fig. [TJ corresponds to the initial conditions for the 
BK evolution (R S (Y) / R s (Yq) = 1), while the right panel corresponds to a higher rapidity Y > Yq 
where R s (Y)/R s (Yq) = 0.34, i.e., after running the evolution for some time. 

For fixed parent dipole size \x — y \ the most extreme situations arise when we vary z' along the x — y 
axis and the axis perpendicular to it, all other directions in the z'-plane interpolate smoothly. The 
two axes, along with one intermediate 45° axis, are also shown in both panels of Fig. [TJ Evolution 
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Fig. 1: Contour plot of the correlators on the right-hand side of the BK equation |?J) (divided by 
(S X y) to normalize the large z asymptotics to —1) in the z' -plane at different stages in the evolution, 
for fixed \x — y\ = R s (Yq). Here z' = (z[,z 2 ). The left panel shows the initial conditions for the 
evolution, while the right panel displays the evolved distribution (see text). The dots mark x and 
y, the locations of the parent q and q. These points always fall on the boundaries between positive 
and negative contributions to the right-hand side of ([7]) ( marked by contour lines going through the 
dots). The half-rays denote the angles at which the factorization violations will be plotted in Fig.\^ 



speed, in cither JIMWLK or BK, is then a consequence of a numerically delicate balance of the 
negative and positive regions of the quantity plotted in Fig. [TJ Since we are talking of evolution 
for S in which generically S is driven to smaller values at fixed dipole sizes as rapidity Y increases, 
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the negative regions in Fig. [T] push evolution forward (these contributions are generically those at 
large \z'\), while the positive regions in Fig. Q] (generically near x and y) slow it down. At fixed 
coupling, any change of evolution speed can be mapped onto a change of relative weight of these 
two contributions. Starting from a non-scaling initial condition like the GB-W model, evolution 
typically speeds up until scaling is reached and evolution speed is maximal. (Scaling here is defined 
as the situation in which all rapidity dependence is carried by the saturation scale so that observables 
like (S r ){Y) become functions of scaling ratios like r/R s (Y) only [60].) This is mirrored perfectly 
in a shrinking of the positive regions from the radially symmetric situation of the GB-W initial 
condition (Fig. [TJ left) to a situation in which there are two separate positive regions near the q 
and q positions (Fig.[TJ right). 

The z' plane symmetries of the dipole evolution equations translate directly into the factorization 
violations A xzy (Y) from Eqs. ijHJ, and also the zeroes at the q and q positions carry over. In [29] two 
of us (Rummukainen and Weigert) had observed numerically that all factorization violations tested 
were positive (i.e., qualitatively acted to slow down evolution compared to BK), and unexpectedly 
small, at least in the regions that contribute to evolution: instead of A ~ 1/AT? ~ 10% at N c = 3 
one found contributions roughly another magnitude smaller. Fig. [2] re-illustrates the observation 
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Fig. 2: Factorization violations from JIMWLK evolution (scaled up by N%) plotted against varying 
parent dipole size at fixed \z'\ = 0.4 ■ R S (Y). The angles 0° (left), 45° (middle), and 90° (right) 
are the angles between z' and r and refer to rays in the z' -plane as indicated in Fig\l\for one fixed 
\r\ = \x — y\. Shown are three different rapidities each. One obtains a reduction by a factor of 
10 compared to the natural size A J ~ 1/N% (or A J ^ \). This was observed earlier in [29], 
with any differences being due to the slightly different correlator geometries chosen here for ease of 
comparison with the discussion below. Shown are "typical" regions that contribute to evolution, see 
Sect. for details. Note also that only the 0° ray shows special structure since it contains strict 
coincidence limits (i.e. the limits where x, y or z overlap), all other angles are qualitatively well 
represented by the 90° case. 



of [29] , namely that the factorization violation A is about an order of magnitude smaller than the 
naively expected 10% in regions relevant for evolution. Fig. [2] plots A xzy (Y) (henceforth referred 
to as simply A without the arguments) from the numerical solution of JIMWLK evolution equation 
which we will describe below. A xzy (Y) is plotted in Fig. [2] as a function of the parent dipole size 

r = x — y (10) 
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for fixed \z'\ with the angle between z' and r being 0° (left panel), 45° (middle panel) and 90° 
(right panel). These directions were also shown in Fig. [I] In each panel of Fig. [2] the factorization 
violation A is plotted for three different rapidities: Y%, Yx + 3, and Y\ + 5, with the exact numerical 
value of Y\ being irrelevant here (along with the value of the fixed coupling constant a s ), as our 
goal in this Section is only to demonstrate the size of the typical factorization violations. 

While the correlator geometries in Fig. [2] differ slightly from those shown in [29] , the magnitudes 
are comparable. We see again that instead of naively expected AiV| ~ 1 one gets AN% ~ 0.1, 
which is an order of magnitude smaller. As rapidity increases beyond the values shown in Fig. [21 
the factorization violation A does not grow significantly beyond the values achieved in the figure. 
In [29] no attempt was made to clarify in which regions of configuration space the factorization 
violation A is small, and no generic discussion of relative importance of configuration space regions 
was given. To fully understand the statement of the smallness of corrections one must expand on the 
discussion given there and first gain a better understanding of where to expect sizable contributions, 
since mapping out all of the configuration space in z', x — y and Y is otherwise not feasible. This 
will also provide the underlying reason for the observed smallness. 



3 Origin and smallness of the factorization violation: an in- 
terplay of saturation and coincidence limits 

Smallness of the specific factorization violations JSJ) are only one facet of a more generic question: 
what kind of deviations from full JIMWLK evolution are caused by the factorization assumption j6]) 
with its associated truncation of the Balitsky hierarchy of the quark dipolc operator? 

Full JIMWLK evolution does not only couple in a full hierarchy of evolution equations for the quark 
dipole operator, it has an even wider scope: It consistently incorporates hierarchies based on any 
n-point correlator. Examples for such distinct hierarchies are obtained by considering the infinite 
set of dipole correlators labeled by all finite dimensional unitary representations 1Z. Each of them 
has its own distinct evolution equation, that can be summarily written as 

— (t r (U x U y ))(Y) = -± Jd 2 zlC xzy [{[U z ] ab tr(t U x t U y ))(Y)-C n (tr(U x U y ))(Y)j . 

(11) 

Here U refers to the group element in the representation 1Z, with analogous notations for the trace, 
generators and conjugate representation. C-ji denotes the second Casimir of the representation of 

the dipole, i.e., for the qq correlator of the BK case it equals Cf = N 9 C N 1 or for a gg dipole it would 

be Ca = N c . The gluon produced in the evolution step is denoted by U z and is of course always in 
the adjoint representation. 

The hierarchies based on the dipole equations (|11|) are by no means all independent (group con- 
straints and coincidence limits may reveal that the Y dependence of the same multi-t/ correlator 
does appear in several hierarchies), nor do they exhaust all the information contained in JIMWLK 
evolution (for instance operators with non-vanishing triality arc absent from the family of dipole 
hierarchies). What is important here is that JIMWLK evolution treats this multitude of correlator 
equations consistently - as long as no truncation assumptions are made. 
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The BK approximation greatly simplifies this intricately interlinked set of hierarchies and may not 
capture all of its features in the process: since there is no simple generalization of the Fierz iden- 

tity ([3]) for ( [U z \ tr(t U x t U y )) (Y) in an arbitrary representation 1Z, it may become impossible 
to consistently generalize the BK approximation to even this class of equations, despite the fact 
that one can write expressions for the BK (large N c ) limit of generic dipole operators. (For the 
cases of 1Z being fundamental or adjoint representations the BK approximation is indeed possible 
and is done routinely.) To consistently include all equations (TTTj) is to go at least one step beyond 
the BK approximation and below, in Sect. 21 we shall see that this can indeed be achieved quite 
elegantly. 

The key feature satisfied by JIMWLK evolution that is violated by BK factorization beyond the 

leading- iV c limit is a set of group identities for the three point correlators ( [U z \ tr(t U x t U )) (Y) 
on the right-hand side of the evolution equations (TTTj) . In what follows we will generically use the 
term coincidence limits to refer to the limits were any pair of points (x and y, x and z, or z and y) 
or all three of them (a;, y, and z) coincide with each other. At the coincidence limits the correlator 

\[U Z \ tr(t U x t U y ))(Y) should inherit relationships that JIMWLK evolution respects on the 
operator level (see Appendix |A1 for their derivation): 

lim [U z ] ab tr(i U x t U y ) = Cnj^tr [U z Ui) , (12a) 

n _ n „J. H Tt a Ti ■Rpn't K K ret 

lim [U z ] tr(t U x t U y ) =C n tr(U x U y ) , (12b) 

z — >y or x ^ 

lim [UJ tr(t U x t U y ) =C n d n , (12c) 

where diz stands for the dimension of the representation (df — N c for the fundamental representa- 
tion, d,A = N% — 1 for adjoint, etc.) and tr denotes the trace in the adjoint representation. While 
the third statement is merely a normalization statement, the first two are remarkable: we read off 
that in the limit of small parent dipole the three point operator on the left-hand side reduces to a 
gluon dipole, no matter what representation 1Z refers to, while in the limit z — > x or y it reduces 
to an 7?.7^-dipole. The latter, (|12b[) . is crucial to ensure the real virtual cancellations in (TTTj) . 



For correlators, the implications of (|12[) go far beyond the isolated points featuring in the limits 
shown. Since the correlation (saturation) length R s = 1/Q S is the only dimensionful parameter, 

_ a b n iz a -R. -regret 

the only scale in the problem, (TT2"|) determines the generic behavior of {[U z \ tr(i U x t U y ))(Y) 
in all of configuration space. 

Configuration space is first divided into two classes in which r is either smaller or larger than 
R s . For each of these classes one has to distinguish two cases according to whether the distance 
between the gluon and the nearest quark is larger or smaller than R s . The configurations are shown 
in Fig. [3] and exhaust all physically different situations (labels "a" through "d" in the figure are in 
correspondence to the equation labels in (fT3"|) below). One infers from (TT2")) that in regions "a" and 
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\r\ < \z'\ 



region "a" 



r\,\z'\<^R s 
region "d" 



H » R s , \z'\ « 
region "b" 



\r\/2 




\r\ >i? s , |z'| *j6 |r|/2 
region "c" 



Fig. 3: Schematic representation of physically distinct cases for correlator configurations. The 
shaded areas denote regions of size R s , the particle coordinates x, z' , y are at the corners of the 
triangles with the gluon position z' marked by a dot. 



"b" the falloff is dipole-like 
|r| <C Rs < \z'\ : 
\r\ «2|z'| >i?. s : 



([U z ] ab tr(f tf.i tf w ))(Y) «C K ^(tr([/ z [/t))(F) , (13a) 

r _ - „/, W 7i° It -rP Tit K K tt 

<[[/«] tr(f [/ x t C/ y ))(r) «CR(tr(^^))(y) , (13b) 



i.e., {[U z \ tr(f U x t U y ))(Y) vanishes like a gluon or 7?.7?.-dipole where \z'\ ^> R s ^> \r\ or 
|r| ss \z'\/2 3> i? s (in the latter region the gluon is near either the q or q, implying either R s ^> 
\z — y\ or \z — x\ respectively. This only occurs when the angle between r and z 1 is near 0°). 
It also vanishes trivially in region "c", where |r| 3> R s , and the gluon is far from both q and 
q (\z'\ pfi \r\/2). In this region all three inter-particle distances are large and force exponential 
suppression, although (|12p gives no additional information about the falloff, leaving us with 



\r\ »i? s .,|z'| £ \r\/2 



([U z ] ab tr(t U x t U y ))(Y)^0. (13c) 



This leaves only one region, labeled "d" in Fig. [31 in which the contributions are not suppressed. In 
this remaining region, all scales are small simultaneously, as one would naively expect in a system 
with a finite correlation length R s . In region "d" we have 



\r\<R, 



([U z ] ab tr(t U x t U y ))(Y)<C n d n . (13d) 



Fig.Q]illustrates this theoretical discussion with contour plots of the three point function as obtained 
from actual JIMWLK simulations. The plots show dependence on qq separation \r\ and the distance 
of the gluon location with respect to the qq midpoint \z'\, with z' perpendicular to and parallel 
to r, i.e., along two of the lines indicated in Fig. [TJ One may notice that the contributions on the 
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axes, \r\ = and \z'\ = have no angular dependence: the first corresponds to zero size parent 
dipoles in which case r does not single out any direction to refer to, the second keeps z firmly in 
the middle of the qq pair while varying its size so that again the angle does not play a role. 



{[U z ] a \r{t a Ujul))/C s d s 
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0.49 - ^ 
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Z = 0° 



([U z '] a "tr(t a U x t b Ul))/C f d f 




0.91 



0.70 — 



0.49 — 



0.28 



0.07 



Z = 90° 



Fig. 4: Behavior of three point correlator as discussed in Eq. (|13[) taken from our numerical solution 
of the JIMWLK equation at some intermediate Y . z' is varied along rays of fixed angle with respect 
to r (0° and 90° ), c.f. Fig. QJ The regions are labeled in correspondence to Fig. \3i the correlators 
display the generic behavior anticipated in Eq. |T5J) . Note that region "b" is only present near 0° 
and completely disappears for 90°. 



Eqs. (|12p and (T3"|) represent but one example of a much larger set of group constraints for more 
complicated correlators that are all inherently true in the full JIMWLK setting, but broken by the 
BK factorization. The BK truncation is geared towards quark dipoles (where 1Z is the fundamental 
representation) , where it approximates the Fierz identity (J3j) by dropping the 1 /N c term 



[U z ] ab 2tv(t a U x t b Ul) « tr(U x Ul)tr(U z Ul) + 0(1/N C ) 



(14) 



This distorts the coincidence limit of the quark dipole version of Eqs. (fT2| into their large N c 
approximations 



hm tr(U x UlMU z Ul) = |tr([/ x [/]) | 2 
tr(U x Ul)tv(U z Ul) = tv(U x Ul)N c 



lira 

->7/ or x 



lim tv(U x Ul)tr(U z Ul) = iV r 2 , 



*x; y — >x 



(15a) 
(15b) 
(15c) 



i.e., it approximates the gluon dipole operator on the right-hand side of Eq. (|12a[) by the square of 
the quark dipole operator, and replaces the constants in the remaining equations by their large N c 
counterparts. 

For correlators, the implications of (fT5|) mirror the conclusions drawn in Eq. (fT"3|) up to corrections 
of order 1/-/V 2 , hence one naively expects the factorization violations to be of precisely that order, 
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unless there is a stronger cancellation at work in the coefficient of that 1 /N 2 term. Eq. (flU)) contains 
all that is needed to assess this issue if one uses ([3]) to recast A in terms of two and three point 
correlators only (the first two terms represent the unfactored correlator): 



±xzy\ 



CH = ^2 ([U z ] 2tr(t a U x t b Ul)) + 



<tr(E4E/t)} 



N r . 



(tr(U x Ut))(tT(U z Ul)) (Y) . (16) 



The four distinct regions of Fig. [3] and Eq. (|13[) can then be addressed in turn (all individual 
correlators are real and positive): 



Region 



\z'\ 3> R s S> \r\: Both the first and the last term inside the brackets of (fTo| are 



exponentially small, but the second term approaches 1. In the extreme case \z'\ — ► oo, \r\ — > 
one finds A xzy (Y) -> \jN 2 c . Contrary to ([U z ] ab 2 tr{t a U x t b W y )), {tx{U x Ul)tx{U z U y )) con- 
tains a ^-independent additive term that survives this limit. If region "a" were to contribute 
to evolution at all, this would destroy infrared safety of JIMWLK evolution (see below). 

• Region "b", \r\ m2\z'\ 3> R s (gluon near q or q): Since there are always two large distances 
involved, all three of the terms in (|16p are exponentially suppressed and the contribution is 
naturally much smaller than l/N 2 . 

• Region "c", \r\ 3> R s , \z'\ \r\/2: With all three inter-particle distances large, all terms are 
exponentially suppressed individually, rendering their sum much smaller than 1/N 2 . 

• Region "d", \z'\, \r\ <i? s : The terms inside the brackets are order N% — 1, 1, and AT? re- 
spectively. Moreover, in the strict coincidence limit x = y = z, they cancel exactly! This 
guarantees a very strong (albeit not exponential) reduction of the coefficient in front of 

The cancellation is slightly less pronounced farther from exact coincidence, for scales |z'|, |r| 
of order R s , before large distance damping at the boundary to the previous regions sets in. 

One concludes that A is strictly bounded from above by 1/N%, but there is only one region left in 
which this bound is actually reached - region "a" . In all other regions strong cancellations reduce 
the contributions to values significantly below this bound. 



So far we have used general arguments based on coincidence limits (| X2[) and on the effects of 
saturation on the dipole scattering amplitude (two- Wilson line correlators) to argue that A is in 
fact much smaller than 1/N% for much of its phase space. To understand the impact of A on the 
evolution let us rewrite ([4]) using |8a|) 



dY 



0*v)OO 



u s N c 
2tt 2 



d 2 z /C 



xzy 



(S XZ )(Y) (S zy )(Y) - (S xy )(Y) + A xzy (Y) 



(17) 



Eq. (|17p shows how A enters the full, untruncated evolution equations. 

As we saw above, somewhat surprisingly, region "a" where the maximal possible factorization 
violation occurs is characterized by small parent dipole size \r\Q s — > but with the gluon produced 
far away, with [z'lQg 3> 1. This region, however, has no impact on evolution at all: it is completely 
power-suppressed by the evolution kernel ([5]) in (TT7|) which goes to zero as the sixth power of 
distances involved: 



(x - yf 



'■(z-y) 



(z'y 



o . 



(18) 
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Were it not for this kernel suppression, JIMWLK evolution (represented via the qq Balitsky hier- 
archy) would receive large distance contributions from this region - infrared safety would be lost. 
Notably, region "a" is the only large distance region that requires suppression from the kernel. The 
remaining large distance regions ( "b" and "c" ) are exponentially suppressed on the correlator level 
and automatically decouple from evolution. That leaves the last region (region "d" ) with its strong 
cancellation of contributions as dictated by the properties of the coincidence limits: it remains as 
the sole channel through which the non-factorizcd contributions affect the energy dependence of 
the qq-dipole. It is this region which was quoted in [29] to contribute the "typical" factorization 
violations without connecting this to the coincidence limits (fT2")) . 

While the argument given in this section does not give a parametric estimate for the size of the 
factorization violations, it does explain why the contributions are naturally much smaller than 
We see that factorization violation A xzy is bounded by 1/N% from above. However this 
1/N% value is reached only in a small subset of (x,y,z) configuration space (in region "a"), which 
is suppressed by the evolution kernel. The relative suppression of the integral of the factorization 
violation over all z's in (JT7J) compared to the first term on the right-hand side of the equation is 
therefore much stronger than the 1/JV| one would naively expect. Note that the generic arguments 
given here also do not allow to determine the sign of the contribution and thus do not allow to infer 
if one should expect JIMWLK evolution to be slower or faster than the factorized BK truncation. 

We might now just push ahead and map out configuration space of the JIMWLK 3-point correlators 
using numerical results from our simulations, to systematically supplement the numerical results of 
[29] and Fig. [2] with contribut ions from the regions not shown there (numerical results will be shown 
in Figs. E] and [6|). Let us instead first give a simple generalization of BK factorization that treats the 
set of equations (jTTJ) consistently and respects the coincidence limits (fT2j) . This generalization will 
restore at least part of the true factorization violation and respect the configuration space pattern 
deduced from (|13[) and (fTB|) . This will likely improve agreement with JIMWLK evolution and give 
some insight into the question in which direction evolution speed is changed by the factorization 
violations. 

4 Gaussian truncation of JIMWLK 

4.1 A step beyond BK 

Our argument for the suppression of 1 /N£ corrections in the previous section were based on satu- 
ration effects and coincidence limits. We observed that the BK equation, while incorporating the 
saturation effects, violates the coincidence limits at the subleading N c level. 

Approaches that both incorporate saturation physics and respect the coincidence limits of the 
general argument given in Sec.[3]are well established in the literature. They take the form of variants 
of the McLerran-Venugopalan model [7-9] and the closely related Glauber-Mueller approximation to 
high energy scattering [2-4] , which can be rigorously established by summing QCD diagrams without 
taking into account small x evolution. All these descriptions fall into a class of approximations of 
the JIMWLK average (. . .)(Y) over Wilson lines U x , that is characterized by a longitudinally local 
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Gaussian averaging procedure that can be cast as 

Y 

(. . ->(n - e*p{-lJ<irjM, Gr.„ ^ }... . (19) 

We refer the reader to [27] for a discussion of how various well known models can be recovered from 
the generic form shown in Eq. (|19[) by choosing specific expressions for Gy> ,xy In the quasi-classical 
limit G encodes a two gluon exchange with the target in the f-channel. That this same generic 
approach automatically satisfies the coincidence limits has also been demonstrated in [27] for the 
case of 7t being the fundamental representation. 

Stepping beyond the quasi-classical limit in [43], Kovner and Wiedemann have suggested an all 
N c evolution equation that merges BK principles with the Gaussian treatment of correlators incor- 
porated in Eq. (|19p by what amounts to applying the averaging prescription to quark and gluon 
dipoles. 

In fact, Eq. Q19p allows us to extend the treatment of [43] beyond the specific case of quark (funda- 
mental) and gluon dipole evolution. This results in a self-consistent treatment of the evolution of 
all generic dipole operators in which one replaces the q and q by colored objects in arbitrary repre- 
sentations 1Z and 7t. Doing so, one finds completely generic expressions for the previously discussed 
correlators (see Appendix IB1 for calculational details and also [12,40-43] for similar calculations): 

(tr(U x U y ))(Y) = d n e - c ^*y (20a) 



([U z ] ab tr(TuJ U y ))(Y) = C n d n e-^Wv.-.+Cv.-B-ffy.-.J-OROr,., ( 2 0b) 
For convenience we have introduced 

Y 

Qy,x V '■= j dY' (Gy' ,xy — ~(Gy',xx + GY',yy)^J (21) 

to denote the combination in which the t-channel gluons enter these expressions. Note that Qy.xx = 
as required by consistency in (|20a.|) . Quick inspection reveals that (|20b[) indeed complies with (|12|) 
as advertised. In fact, this property is not specific to this particular set of correlators. Any correlator 
calculated using (flU)) (or any generalization thereof) will automatically satisfy all necessary group 
constraints by construction. 

This procedure then is a candidate to generalize the BK factorization in which one simply trades 
an evolution equation for (S xy )(Y) for an evolution equation for QY,xy The procedure at least 
qualitatively repairs the flaw that is the source of factorization violation in the BK equation. We 
will find below that it provides quite good qualitative insights on factorization violation but is not 
sufficient to obtain quantitatively correct results. 

The equation for Q has already been derived in [27] (and in a somewhat different form earlier 
in [43] 5 ), starting from the gg-dipole evolution equation (|2|). We reiterate that this average treats 

5 To connect with the form given in [43], Eq. (5.3), one should reconstruct the evolution equation for the qg-dipole 
operator by multiplying ( 1221 ) with exp(— CjQy.xy ) an d note that our Qy,xy corresponds to v(x, y) in [43]. 
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all dipole equations consistently: Inserting (|20p into the generic dipole evolution equation pip 
yields one and the same equation for Q, 

^Gy, xv = g \d 2 z K xzy (l e -H^*^v-S^ , (22) 

irrespective of the representation 1Z. Note that ([22]) is similar to (though not exactly the same as) 
the Ayala-Gay Ducati-Levin (AGL) evolution equation [61-63]. 

Contrary to BK evolution which systematically discards all 1/N% suppressed terms contained in 
JIMWLK, the Gaussian truncation has no expansion parameter justifying the approximation. Nev- 
ertheless we expect it to lead to a good approximation of JIMWLK evolution since 

• the equation incorporates a subset of these 1/N% corrections that is sufficient to restore the 
coincidence limits; 

• the low density limit of Eq. (|22|) (viewed as its small Q limit) reduces to the BFKL equation; 

• it has a large N c limit that is compatible with the BK equation as will be seen below. 

Despite (|22[) being surprisingly more generic than the BK equation, in the sense that the procedure 
addresses arbitrary dipolcs irrespective of representation, one remains with an approximation to 
the true JIMWLK evolution, and does not obtain an exact solution of the JIMWLK equation: the 

~ a b ft iz a ft ft b ft t 

evolution equation for ( [U z \ tr(t U x t U y )){Y) resulting from JIMWLK would impose additional 
conflicting conditions on Q, and can only be satisfied by introducing degrees of freedom beyond Q. 
The Gaussian truncation still deviates from JIMWLK evolution at the level of evolution equations 
for three point functions. 

It is worth noting two particular features that the average entails. First, Eq. (|20a[) implies 
Casimir scaling for dipole correlators: 6 Given two representations 7?-i ; 2 the normalized dipole cor- 
relators are related by a simple power law 

1 fti fti 7*i t ( \ ft 2 ft 2 «at \ C ^i/Cft 2 

— (tr(U x U y ))(Y) = ^—(tr(U x U y ))(Y)j . (23) 

Second, somewhat surprisingly, Eq. (|22[) can be mapped back onto the BK equation. This implies 
that the dynamical content of the Gaussian truncation is the same as that of the BK equation. The 
main improvement is how this information is mapped onto the correlators. As we shall see, this 
leads to a slightly better approximation of JIMWLK results. On the practical side, this turns into 
a time saver: one can recycle the numerical tools written to solve the BK equation, provided one 
relates correlators and initial conditions accordingly. 

One way to see that the dynamical content is the same is based on a simple re-parametrization of 
the BK S-matrix in as close an analogy to (|20a[) as possible. We write 

1 - N™ y = = MU x Ul)){Y)/N c = e -*^.-» . (24) 



6 Approximato Casimir scaling has been observed for Wilson line correlators in the context of heavy quark poten- 
tials in [64,65]. 
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Qy,xy should be thought of as simply being defined by the solutions of the BK equation through (|24D . 
The N c dependent constant is a convention chosen in keeping with the N c lore of the Mueller dipolc 
model and the BK equation. Next one inserts this into the BK equation for S, 

J_cBK _ a sN c r 2 / BK BK BK \ . » 

Y > X V ~ 2lT 2 ^xzy y^Y,xz D Y,zy s Y,xy J l \ z0 ) 

and obtains 

^Gy,x V = % J*z K xzy (l - < <^ - (26) 

which is identical to (|22|) . thus establishing our claim of identical dynamical content, despite the 
different treatment of correlators. 

The procedure to obtain solutions for the Gaussian truncation that would serve to determine, say, 
the gg-proton cross section for DIS at HERA would then be to choose an initial condition for S, 
read off Q via (|20a[) . insert it in place of Q in to determine the initial condition on S' BK to be 
used in the BK equation (|25|) . After solving Eq. (|25|) to obtain S BK at all rapidities one reverses the 
procedure to recover Q via ([24]) and (|20a[) at each rapidity Y. This Q then determines all correlators 
of the Gaussian truncation through (fT9|) and the special cases shown in Eq. (1201) . 

As an immediate consequence we conclude that the asymptotic scaling shape for the dipole correla- 
tors in the Gaussian truncation can be obtained from those of the BK equation using a simple power 
law relationship l|23p . As an immediate corollary also evolution speeds of the Gaussian truncation 
and BK evolution must coincide in that region. This link does not extend to the pre-asymptotic 
regime and we will see that GT tends to be slower that BK in that range in Sec. O 

Two further points are worth noting: (i) one may recover BK factorization (wherever it can be 
meaningfully applied) as the leading N c contribution of the new procedure. (This can be verified for 
the contributions entering the BK equation by taking the large N c contributions in the exponents 
in (|2"0)) .) (n) In the small density limit, i.e., the limit of weak target fields where Q is small its 
evolution equation, Eq. I|22[) . consistently reduces to the BFKL equation for Q, in keeping with the 
underlying interpretation. 

In summary, one might think of the Gaussian truncation as a minimal extension of the BK factoriza- 
tion to consistently incorporate group constraints with an associated set of "minimal" sublcading 
1/N C corrections without changing the dynamical content of the associated evolution equation. 
Below we will frequently abbreviate Gaussian truncation as GT. 

4.2 Factorization violations in the Gaussian truncation 

Now that, with JIMWLK, BK and GT, we have accumulated three different ways to simulate small 
x evolution for all of which the averaging procedure (...)(Y) leads to different results we need 
to refine our notations to avoid confusion by distinguishing (...)j(Y), (...)b(Y) and (. . .)q(Y) 
respectively. We also are faced with different factorization violations and define 

A izy( Y ) ■■= {{Sxz - (Sxz).i(Y)){S zy - (S.yhOO))^) (27) 

to distinguish it from the analogous (non-vanishing) quantity taken in the Gaussian truncation 
which we will denote A£ (Y). Only in BK this is set to zero by fiat: (Y) = 0. 
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A£ zy (Y) is the only channel through which higher order correlations of JIMWLK feed into the 
dipolc equation of the Balitsky hierarchy as can be made explicit by splitting ((4]) into two parts as 

J^(S xy ).j(Y) = J d 2 zK xzy ([S xz )j{S zy )j - (S X y)j^j (Y) 

+ ^ Jd 2 zK xzy Ai zy (Y). (28) 

Without the A' 7 term on the right-hand side this would reduce to the BK equation and decouple 
from the rest of its Balitsky hierarchy. As such it is also the only source of possible differences in 
evolution speed and correlator shape between JIMWLK, BK and the Gaussian approximation. The 
evolution equation of the latter can also be rendered in the form (|28|) with all J replaced by G. 

The Gaussian truncation provides explicit expressions for its associated factorization violation 
A^ zy (Y) which manifestly respects the regional patterns outlined in Sec. [3] on general grounds. 

Using (f2"0|) and the Fierz identity © one finds 

A G y (Y) = -!d J { —£ il e- c '<*- , (29) 

where By,xzy is a shorthand notation for the correlator combination already known from the evo- 
lution equation ([22]) . the expression for the three point function (|20b[) or the right-hand side of the 
BFKL equation 

By,xzy '■= Qy,xz + Qy.zy — Gy,x V • (30) 

A G is positive semidefinite and strictly vanishes only where B vanishes. B = is the hyper-surface 
on which the integrand of (|22[) changes sign, in analogy to what was shown for the BK equation 
in Fig. [TJ At large positive B, the fractional expression in (|2"9")) approaches 1/N%, and it grows 
exponentially at negative B. The region of large, positive B corresponds to \z'\ \r\ (region "a" in 
the notation of Sect. [3j, the region power suppressed by the evolution kernel. On the other hand, 
the region of maximally negative B corresponds to large parent dipolcs |r| 3> R s (regions "b" and 
"c" in Fig. [3]) and is is strongly suppressed by the overall e~ CfSY ^ y . This leaves the region near 
B = to contribute. This region coincides with the region "d" singled out in Sect. [3] to yield the 
main contributions. Hence A G in (|29[) illustrates the main feature of the factorization violations 
derived in Sect. [3]). 



The factorization violation in the Gaussian truncation Eq. (|29|) also explicitly vanishes in the low 
density (BFKL) limit: Contributions to Eq. (|29|) in fact start off at order B 2 , i.e., arc manifestly a 
nonlinear effect. 

Since the right-hand side of the BK equation for (S) (0 is negative, one immediately concludes 
that the positive contribution of the factorization violations included in the Gaussian truncation 
slow down evolution for any parent dipole size r. One would expect this to carry over to JIMWLK 
evolution, i.e. one would expect JIMWLK evolution to be slower than BK. 

In the following we will show plots that map out factorization violations and their influence on 
evolution. Since we want to show both the dependence on parent dipole size \x — y\ and the position 
of the produced gluon z we can not simultaneously plot A against all z degrees of freedom. We 
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rely on the z'-plane symmetries shown in Fig. [T] and restrict ourselves to half rays at 0° and 90° 
as done earlier and remind the reader that other half rays will smoothly interpolate these extreme 
cases. 

Contour plots of the factorization violations in the Gaussian truncation and the full JIMWLK 
simulation in the continuum limit (see Sect. [5] below) at some fixed rapidity are shown in Fig. [5j 
They confirm that the regional pattern deduced in Sect. [3] is indeed seen in both JIMWLK evolution 
and the Gaussian truncation. Despite this qualitative agreement, it becomes evident that the 
Gaussian truncation underestimates the magnitude of the contributions significantly. Fig. [6] shows 
the corresponding contributions to the integrand of the evolution equation (|28j) . i.e. after the kernel 
has been multiplied in. The kernel both power suppresses the large \z'\ region and enhances the 
contributions of small \z'\ and, as already anticipated, the results become quite sensitive to the 
short range behavior as B goes to zero in the lower left hand corners of the plots. 



5 Quantitative consequences: slowdown of evolution 

The notion that the presence (or absence) of a factorization violation term in Eq. (f2"5| would affect 
evolution speed can be made more precise. By extension of an argument in [66] we define evolution 
speed 7 as a J d 2 r/r 2 -integral of the right-hand side of the evolution equation 

A(F) := /S" Kxzy {§xz§zy ~ §xy)iY) (31) 

and then proceed to split the contributions according to (|28|) . (Again r = x — y.) 

To obtain a quantitative comparison of JIMWLK and BK equations, one should take note that due 
to non- vanishing A J 

(S xy ).,(Y) ± (S xy ) B (Y) (32) 

even if one chooses them to be equal at the initial rapidity Yq\ 

(S xy )j(Y ) = (S xy ) B (Y ) (33) 
as will be the case in all the numerical comparisons shown. 

To calculate the difference of evolution speeds at some finite Y — Yq one has to compare Eq. (|28[> 
with the right-hand side of 

^(S XV ) B (Y) = ^r/ d 2 zlC xzy {(S xz )b(S zv )b (S xv )b) (Y). (34) 

The difference between the JIMWLK and BK evolution speeds from Eqs. (|28|) and (|34[) is due to 
the difference in the JIMWLK and BK averagings of the S-matrices and to the presence of the 

7 Evolution speed as defined in [66] refers to the scaling regime with a uniquely defined Q B (Y) and is defined as 
A(Y) := In Q 2 S {Y). Outside the scaling region the starting point Qs{Y) is no longer uniquely defined and one 
may take this as one of many possible definitions for a well behaved measure of evolution speed. 
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Fig. 5: Contour plots of the factorization violation A scaled up by at a fixed rapidity in 
JIMWLK (left) and in the Gaussian truncation (right). The plots scan parent dipole size \x — y\ 
and distance \z'\ from the midpoint (x + y)/2 at 90° (top row) and 0° (bottom row) with respect 
to x — y. The 0° case is special since it contains contributions where A strictly vanishes due to 
the coincidence limit z = x or y which appears here as a line with \z'\ = \x — J/|/2. The bulk of 
the contributions is similar to the 90° case. As was the case for the three-point-functions of Fig. ^ 
the contributions along the axes \r\ — and \z'\ = are identical for all angles. Fig. [H shows cuts 
along horizontal lines near the bottom of the two JIMWLK plots. 
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Fig. 6: Factorization violations shown in Fig. \3\ multiplied by the BK kernel (i.e. the integrand of 
the A J term in the evolution equation S28\) that determines how much they contribute to evolution), 
by Q 2 S N 2 , and by a scaling factor with varies from panel to panel. The notation and geometry are 
the same as explained in the caption of Fig. [31 In both cases the region where A reaches 1/iV^ is 
completely suppressed by the kernel. Note the marked difference in size of the rescaling factor: the 
contribution of factorization violations in JIMWLK (left, with 10 4 J is an order of magnitude larger 
than in the Gaussian truncation (right, with 10 5 ,). 
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A J -term in (f2"5)l . Indeed these two causes of difference are interconnected. Introducing the shape 
dependent correlator difference 

A£* (T) := ((S xz )j(S zy )j - (S xz ) B (S zy ) B - (iS xy )j - {S xy ) B )) (Y) (35) 

the difference in evolution speed arises from two non-vanishing contributions according to 

A\ JB (Y)~\j(Y)-\ B (Y) (36) 

= - ^/pF / ****** *UY) ~^J^J ***** ^) ■ 

Given identical initial conditions, (|33|) . the A J -term is the sole reason for a non- vanishing AA ,/B (Y) 
to be generated at all, but the second term may in some cases be quantitatively not less important 
once that has happened. 

It is worth noting that in the definition of A, Eq. (f3Tj) . the contributions at small r 2 are enhanced 
by the conformal measure d 2 r/r 2 . This carries over to the differences of evolution speed between 
factorized and unfactorized evolution, Eq. (|36|) . Wc have already shown in Fig. [6] that this is where 
the main contributions to the product of kernel /C times A lie, so that one might expect that AA 
receives an enhanced contribution from any non- vanishing A observed in our simulations. 

Since it is the product of kernel and A that enters one would expect that A should also be sensitive to 
any modifications to the kernel as wc step beyond leading order by, for example, inclusion of running 
coupling effects. We would expect running coupling effects to gencrically suppress contributions at 
small \r\ and thus reduce some of the speed difference visible at leading order. At the same time 
one should be aware that contributions beyond running coupling would also affect evolution speed 
and the way 1/N C corrections couple into dipole evolution. Those might be of the same size as 
the running coupling contributions, but it would be very peculiar if they were to affect the general 
pattern of the leading oder behavior observed in the simulations shown below. 

We have performed a numerical solution of the JIMWLK evolution equation at one loop accuracy 
with fixed coupling along the lines of and using the techniques developed in [29]. The method is 
based on the fact that the JIMWLK equation takes the form of a functional Fokkcr-Planck equation 
which has an equivalent Langevin formulation [18,67] in which the averages (. . -}j(Y) are expressed 
as ensemble averages of fields Uy, x , with the Y-dcpcndcncc implemented by a Langevin equation 
for the ensemble members. To implement such a description numerically one is forced to discrctize 
transverse space to represent the ensemble member fields Uy, x in terms of a finite number of degrees 
of freedom. This automatically introduces a UV regulator in the guise of the lattice spacing a and 
a IR regulator, the lattice size L. The reader interested in further details of the numerical solution 
and in the implementation of the numerical procedure is referred to [29]. In order to be able to 
directly compare our JIMWLK simulations with BK and GT results at fixed a and L, we have 
chosen to run our BK and GT simulation on the same type of regular lattice with identical a and 
L values although that is decidedly not the most efficient way to perform standalone BK or GT 
simulations. 

Fig. [7] shows numerical results for evolution speeds A(Y) as a function of R S (Y) = 1/Q S (Y), for 
lattices sizes varying from 256 2 to 4096 2 . The JIMWLK results of the first two plots are compiled 
from the simulations presented in [29], this were the largest lattice sizes taken into account then. 
All the others are based on new simulations, with measurements taken at smaller Y intervals (which 
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Fig. 7: Comparing evolution speed in JIMWLK, BK and GT on different size lattices, starting 
from initial conditions with identical dipole correlators {S xy ) j(Y<f) = (S xy ) b(Yq) = {S xy )GO / o). 
The vertical line in the first four plots marks L/R S (Y) = 20. The bands show statistical (Jackknife) 



explains the more fine-grained raggedness of the JIMWLK results - adding additional intermediate 
steps would further enhance the phenomenon). 

For the new runs, initial conditions were chosen to carefully explore the convergence to both the 
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continuum limit and the infinite volume limit, i.e., to scan for UV- and IR-cutoff artifacts. IR phase 
space available in the simulation at the initial condition at Yq is varied by increasing lattice size L 
compared to initial correlation length R s (Yq), UV phase space is varied by increasing lattice size 
at (approximately) fixed L/R s (Yq). During evolution active phase space, which is centered around 
1/R S (Y), moves towards the UV, so that available IR phase space grows with L/R S (Y) while the 
UV shrinks with R s (Y)/a. 

The most striking feature of the JIMWLK simulations are the large fluctuations of evolution speed 
on all but the largest lattices. The fluctuations turn out to be IR dominated, they average out as we 
increase the number of points in the IR. UV cutoff effects manifest themselves only for the longest 
runs, as a relatively sharp turn downwards as one follows the curves from right to left as R s shrinks 
while Y grows. This downturn (where present) indicates that the simulation is running out of UV 
phase space with R s {Y)/a < 10. While this behavior is not physical, it affects all our simulations 
in the same way and one does not prevent us from comparing the behavior of the simulations with 
each other. 

This can be read off from Fig.[7|by tracing the following systematic features: With the saturation 
scale safely more than an order of magnitude smaller than the inverse lattice spacing, only varying 
IR phase space affects A. As L/R s (0) increases from 8.35 to 99.46 JIMWLK evolution becomes less 
and less affected by IR fluctuations. Evolution speeds (as compared to BK and GT, which both 
are not affected by fluctuations) slow down until they settle at their infinite volume limit at around 
L/R s (0) w 50. 

Note that in the first four panels (with smallest L/R s (0)) JIMWLK is initially faster than both BK 
and the Gaussian truncation, before this becomes less and less pronounced as L/R S (Y) grows with 
evolution. This is a direct consequence of the impact of fluctuations becoming less pronounced 
as shrinking R S (Y) cuts off contributions from the infrared. It turns out that for very small 
L/R S (Y) where IR fluctuations contribute most, the A JB contribution completely overwhelms the 
A J contribution which in all cases gives a contribution that slows evolution down. The relative size 
of this contributions shrinks strongly when L/R S (Y) > 20 and is no longer able to overwhelm A J 
in the runs with L/R s {Yq) > 25. 

Similarly, for all other runs, where L/R 8 (Y) > 25 from the outset, we observe (on average) a clear 
hierarchy of evolution speeds with BK the fastest, the Gaussian truncation in the middle and 
JIMWLK the slowest. In this range this is simply a reflection of the relative size of factorization 
violations: the larger these are in the regions enhanced by the kernel, the slower the evolution 
becomes. This hierarchy is already visible in the rescaling factors of Figs. [5] and El scaling violations 
in GT are consistently an order of magnitude smaller than in JIMWLK across the Y range explored. 
This becomes evident again, if we contrast the JIMWLK results of Fig. [2] with their counterpart in 
the Gaussian truncation, Fig. [S] Were we to extend our comparison of BK and GT evolutions into 
the asymptotic range, however, evolution speeds would necessarily become identical as discussed 
in Sec. 14.11 We can use this to assess how closely the simulations shown in Fig. [7] approach the 
asymptotic scaling region. With the regular grid necessary to compare with JIMWLK this is 
not practical, but a simulation that only compares BK and GT can make be implemented more 
efficiently and in fact reach the asymptotic limit. Using this freedom we find that this occurs just 
beyond the region where we loose UV accuracy in the longest JIMWLK simulations such as that 
in the bottom middle plot of Fig. [7] The result of this comparison is shown in Fig. [9l Comparing 
the ratios of evolutions speeds in Fig. [9] we conclude that in the asymptotic scaling region (at fixed 
coupling!) one should expect a factorization violation induced 3-5% slowdown in evolution speed 
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Fig. 8: Factorization violations in the Gaussian truncation (scaled up by N%) against varying 
parent dipole size at fixed \z'\ = 0.4 ■ R S (Y) depicted here for comparison with JIMWLK results 
shown in Fig. [H The larger factorization violations in JIMWLK lead to slower evolution. The 
strong change form Y\ to higher rapidities is mirrored by a convergence of evolution speed between 
GT and BK with evolution towards asymptotic regime shown in Fig. [3 

in JIMWLK evolution compared to BK evolution at one loop accuracy. There is no sign from 
the simulation and no theoretical reason to argue that factorization violations in JIMWLK should 
disappear in the asymptotic region. As already noted we would expect this relative slowdown effect 
to become less pronounced at NLO. 

6 JIMWLK beyond the Gaussian truncation: higher order 
correlators 

The simulation results shown in the previous section show a persistently small but measurable 
improvement of the Gaussian truncation over the BK approximation. Still, JIMWLK evolution is 
much more general than either truncation. Both approximations restrict the information retained 
in evolution to two point functions that, in the low density limit matches up with double reggeon 
exchange as incorporated in BFKL evolution. JIMWLK evolution, on the other hand, allows 
for multi-reggeon exchange in evolution and even the limited set of correlators discussed in the 
above is sensitive to their contributions. An example of this is the simplistic form in which dipole 
correlators of higher representations are mapped back onto the quark dipole correlator. This is a 
direct consequence of the fact that the Gaussian truncation only iterates two reggeon exchange into 
Glauber exponents. 

One might attempt to include multi-reggeon exchanges by generalizing (|19p to include higher order 
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Fig. 9: Ratios of evolution speed with finite lattice spacing and in the continuum limit. The plots 
correspond one to one to the bottom row in Fig. [The sharp up- or downturn of the curves with 
finite a towards the left indicate complete breakdown of the simulations in the UV.] Shown are in all 
plots, from bottom to top (excepting the dash-dotted lines): Xj(Y)/Xb{Y) and Xg(Y)/Xb(Y) (both 
at finite a) and Xg(Y)/Xb{Y) (in the continuum limit). A comparison of the two upper curves 
indicates the size of the UV cutoff effects which remain fairly small compared to the error on the 
JIMWLK results. Evolution in JIMWLK is slower than in both BK and GT in all plots. The left- 
and rightmost plots serve to assess lattice artifacts in JIMWLK evolution. Compared to the middle 
plot, the left plot has the same IR but less UV phase space, the right plot less IR but the same 
UV phase space. Both of these "fall through the lattice" already at R S (Y) / R S (Y ) ps .5, they agree 
with the middle plot within errors. (To facilitate this comparison the solid line from the left panel 
is repeated as a dash-dotted line in the other panels.) Using the middle plot where the JIMWLK 
simulations reach closest to the asymptotic regime (where Xq{Y)/Xb(Y) — ► 1 in the continuum 
limit, top curve), we read off roughly a 3-5% slowdown of JIMWLK- relative to BK-evolution due 
to factorization violations. This is expected to be strongly reduced by running coupling effects. 



terms; naive inclusion of three-reggeon terms would modify (|19[) to 

(...)(F) = exp|-i jdY' Jd*xd 2 yG Y ,, xy -J^-J^- (37) 

and include odderon contributions [49]. However, starting with three reggeon terms, locality in Y 
is an assumption that might prove to be too restrictive and will generally not lead to a consistent 
treatment of the Balitsky hierarchies. Moreover, one has to be careful in simply exponentiating the 
3-rcggeon terms, as is done in (|37j) . Here the best way to keep the calculations under parametric 
control is to employ the power counting developed for the classical gluon fields in [10, 11]. In the 
quasi-classical limit the leading term in the exponent of (|37[) corresponds to a two-gluon exchange 
with a nucleon in the nucleus, such that Gy, X y ~ a ^ A 1 / 3 with A the atomic number of the nucleus. 
For a s <C 1 and 4 » 1 there exists a regime where a 2 s A 1 / 3 ~ 1 and the GT approximation of 
Sect. 2]resums all powers of a 2 A 1 / 3 . From the standpoint of this quasi-classical power counting 
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the second term in the exponent of (|37[) corresponds to a 3-gluon i-channel exchange with a single 
nucleon. such that G Yxyz ~ G Yxyz ~ A 1 / 3 , i.e., it is suppressed by one power of the coupling 
a s compared to the leading term. Iteration of such term more than once would be beyond the 
precision of the approximation: two 3-gluon exchanges arc of the same order in a s and A as a 
2-gluon exchange combined with a 4-gluon exchange. From this perspective the contributions of 
Eq. (|3T[) are only under parametric control up to linear order in G* and G d : 

(. . .xn = cxp { -I ) d r J M y Gy ,„_|__A_ J ,38) 

Y 

1 - j x JdY' J d 2 xd 2 yd 2 z (G Y , xyz f abc + G Y , xyz d abc 



SA X + Y , 5A b + Y , 5A C + Y , 



A discussion of truncations that systematically include multi-reggeon contributions goes beyond 
the scope of this paper, but it is not hard to play with (|37|) as an ansatz to explore the conse- 
quences of the inclusion of an oddcron term in this manner (see Appendix [Cjl: it becomes quite 
manifest that such multi-reggeon contributions naturally break Casimir scaling. The reason for 
this is that in general, higher representation contributions in the i-channel start to pick up on the 
more complicated decomposition of a general s-channel 7?.7?.-dipole into irreducible representations. 
(Higher representations were also considered in [68].) For the odderon contribution they account 
for the fact that the qq dipole acquires an imaginary part (with a specific N c dependence) while a 
gg-dipolc remains real since the adjoint representation is real by definition. 

With the numerical results from JIMWLK at hand it is straightforward look for the actual presence 
of Casimir scaling-violating effects in JIMWLK evolution. This is explicitly shown in Fig. [TDJ We 
note in particular that the violations of Casimir scaling do not grow with energy: they seem to 
qualitatively scale with Q s and might, if anything even be slowly erased, but any firm conclusion 
to that effect is beyond the present numerical accuracy, in particular because of the short lever arm 
available before the simulation starts to "fall through the lattice" (i.e. runs out of UV phase space). 
Note that this qualitative Q s -scaling behavior occurs in a region far outside the Q s -scaling region of 
the dipole correlator itself. Presently we have no systematic explanation for this observation. (The 
Q s -scaling region of the dipole correlator is not reachable in the current simulations of JIMWLK 
equation due to the limited UV phase space on the lattices used.) 

To illustrate that the violations of Casimir scaling are driven by nontrivial coordinate dependence we 
show in Fig. llll that no power law relationship modeled on Eq. (|23|) can provide a good explanation 
for the differences observed in Fig. [TO] By itself, this does not exclude a more intricate functional 
relationship, but that in turn would require its own explanation. Our earlier arguments would lead 
us to believe that it is much more likely and natural that new degrees of freedom (starting with 
four point contributions to the average (|37|) ) are needed to explain this difference. 



7 Conclusions 

We have explored the size and nature of 1/N C corrections in the JIMWLK equation and have found 
that quite natural cancellations lead to the much stronger suppression than the naively expected 
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Fig. 10: Violation of Casimir scaling for fundamental and adjoint two point correlators for three 
different rapidities plotted as a function of dipole size r in units of correlation length R S (Y). Left: 
adjoint correlators compared to the Casimir-rescaled fundamental correlators. Error bands indicat- 
ing numerical uncertainty are too narrow to be clearly visible. Right: differences of the two. (These 
are more stable than ratios which tend to become uncontrollable at large distances.) The bands 
indicate the size of the errors. The violation follows correlation length R S (Y), indicated by vertical 
lines. It would appear to scale with Q s well within errors. Note that this is in a regime where the 
dipole cross section has not yet reached its scaling regime and decidedly does not scale yet. 



1/N% ~ 10% (for N c = 3) observed in [29]. The argument is based first on both group theoretical 
coincidence limits for singlet Wilson line correlators and scaling with the saturation scale Q S {Y) to 
establish suppression of contributions in most of configuration space. All remaining contributions 
are then shown to be then decoupled from evolution after suppression by the BK kernel has been 
taken into account. 

We have shown that 1/N C corrections enter through the factorization violation A and have explored 
its properties. It is bounded by 1/AT? from above, but is much smaller than 1/N% for most of its 
phase space, due to saturation effects controlled by coincidence limits leading to extra suppression. 
The argument is generic and can be easily applied beyond the leading ln(l/x) approximation used 
here. Thus saturation effects provide an extra suppression of 1/N% corrections to the BK evolution 
that reduces the difference between the JIMWLK and BK results for the dipole scattering amplitude 
considerably at any accuracy. While NLO corrections will have some quantitative impact, we do 
not expect them to grossly change the LO result for A or correlator differences (when compared at 
the same R S (Y) or Q S (Y)). At one loop we find typical contributions to A of the order of 10 -3 or 
0.1% of the individual correlator values for correlator differences. 

To complement the qualitative discussion for correlators in JIMWLK evolution not only numerically 
but also analytically, we have made a step beyond the leading 1 /N c BK truncation by exploring an 
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Fig. 11: The violation of Casimir scaling in JIMWLK evolution for two lattice sizes with the same 
infrared phase space as indicated by the L/R S (Y ) value. The plots explore the correlator difference 
^xy ~ {^xy) b as a function of dipole size r for values of b = 2 (the BK result) up to the GT value 
of N c /Cf = 9/4 at N c = 3. Intermediate values of b fall into the shaded areas. The result for "a 
best fit" is indicated as a dashed (red) line. Clearly a simple modification of the power alone is not 
sufficient to remove the mismatch. To remedy the situation one would either need a more complex 
functional relationship between dipole correlators or, more likely, new degrees of freedom to alter 
the (x — y) 2 dependence in at least one of the correlators shown. 



alternative truncation of JIMWLK evolution in the spirit of a Glauber-Mueller iterated two-reggeon 
exchange truncation that we dubbed the Gaussian truncation. This truncation includes a minimal 
set of subleading 1/N C corrections necessary to restore the coincidence limits for correlators that are 
at the core of our cancellation argument for subleading contributions. Correspondingly, it includes 
a minimal set of factorization violations. They turn out to have the right qualitative structure but 
are numerically noticeably smaller than the factorization violation of full JIMWLK evolution. The 
Gaussian truncation allows access to a larger subset of the Balitsky hierarchies than BK truncation 
by treating evolution equations for dipole operators in arbitrary representations consistently. As 
a result the Gaussian truncation proves to be a somewhat better approximation to full JIMWLK 
evolution. Accordingly, one of its main consequences, Casimir scaling between dipole operators of 
different representations, turns out not to be strongly violated in full JIMWLK. 

We have firmly established that factorization violations slow down evolution compared to the BK 
truncation, both in the minimalist from introduced in the Gaussian truncation and in full JIMWLK. 
One of main differences between the two is that factorization violation in JIMWLK persist during 
evolution on a level comparable with the factorization violation present in the GB-W-like initial 
condition used, while they become notably smaller in the Gaussian truncation. 

Evolution speed is somewhat more sensitive to 1/N C corrections than the factorization violations in 
the correlators due to an enhancement of contributions from small parent "dipoles." At one loop 
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accuracy this leads to a relative slowdown of JIMWLK evolution that approaches 3-5% near the 
scaling region. Running coupling corrections are known to reduce evolution speed by suppressing 
the relative importance of small parent dipoles. We have argued that the same mechanism is likely 
to also reduce the relative slowdown, i.e. the impact of 1/N C corrections on evolution speed. 

The subleading-iV c terms present in JIMWLK but absent in BK can likely be attributed to multi- 
reggeon i-channcl exchanges and multi-reggeon splitting vertices. The smallness of the difference 
between the dipole scattering amplitude obtained from JIMWLK and BK appears to indicate that 
multi-reggeon effects are not important for this observable. Further investigation is needed to 
clarify if this is indeed the case. 

Casimir scaling violations present in JIMWLK evolution provide a means to assess multi-reggeon 
exchange contributions. We have illustrated this both with a sketch model that includes oddcron 
contributions and a numerical comparison of gg and qq dipoles. While the oddcron contributions 
play no role for the observables considered here and arc gcncrically suppressed by evolution, the 
multi-reggeon contributions present here are small but contribute throughout evolution. We have 
numerical hints at Q s (Y)-scalmg behavior of these contributions that sets in much earlier than 
geometric scaling of dipole correlators. At present we have no systematic explanation for this 
observation other than "naturalness." 

Our whole discussion was carried out at the leading log level, without any NLO contributions 
taken into account, despite the fact that they are known to strongly influence evolution speed. 
This is partly due to necessity: no numerically practical way has been devised to include higher 
order effects, for example the running coupling corrections obtained in [32-34]. However, since the 
arguments given for the suppression of factorization violations are completely generic, one expects 
the observations made here to persist to higher orders in the perturbative expansion, even though 
quantitative modifications are expected to affect the precise numerical result of the cancellations 
observed at leading order. 
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A Generic coincidence limits 

To understand the coincidence limits for the operator \U Z \ tr(i U x t U y ), generic properties of 
the adjoint representation are useful. They immediately give rise to the identities 

r ~ ,„i. x k"k tzP TZ- 1 _ i tz n nC-Rpn'S _ n iz a -R. ret K fc 

[U t ] tr(t U x t U y ) = [Ulu z ] tv(U x t t U y ) = [U z U y ] tr(t U x U y t ) (Al) 
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by pulling adjoint factors out of the trace. 
In the limit x = y this is proportional to 



tr(t t ) - a n 5 ab , 



(A2) 



and the first task is to understand the constant. Clearly it has to be proportional to the Casimir of 
the representation, but normalization conventions do also play a role: With standard conventions 
the fundamental representation has af un d = olf = ^, while in the adjoint representation one obtains 
^adjoint = a a = N c . This can, in fact, be expressed via the Casimir values and the dimension of 
the representation. The usual definition of the quadratic Casimir, 



1Z a 1Z a 

t t 



implies 



TZ / TZ TZ 



tr It t = C v d 



and thus, together with the definition of an above, and-A = Cndn or 



an = C- 



■r 



(A3) 



(A4) 



(A5) 



with cIa the dimension of the adjoint representation. This readily leads to aA = Ca = N c and 

2N c C f ~ 



ap = Cf — h a s obtained from direct calculation. 



This allows to write the x — y limit as 

,_ ~ _,„), n TZ a TZ -Rpn^ ,__ ~xinr TZ n" TZ TZ t TZ b r ~ _ nr TZ TZ a TZ c 

lim [U s ] ti(t U x t U y ) = lim [U z U y ] tr(i U x U y t ) = [U Z U£\ tr(i t ) 

y >x y >x 

= [u z U x ] ac a n S ac = C K ^tr(u z Ui 



(A6a) 



Note that this leads to a correlator in the adjoint representation, irrespective of TZ. The only 
reference to TZ is the proportionality factor. 

The limit z — > y or x is simpler: 



,._ n „J, TZ TZ a TZ ~rP TZ t TZ TZ a TZ TZ t jz a TZ Tit TC TZ TZ t 

lim [U z ] tv(tU x tU y )= lim tr(t U X U Z t U z U y ) = C n tr{U x U y ) 

z — >y or x u z *y or x y y 



The completely local limit can be obtained directly from (IA6a[) to be Cndn or via 

TZ TZ a TZ. ~Rp TZ t TZ TZ a TZ TZ^ JZ a TZ TZ t 

tr(t [/„) = lim 



(A6b) 



■R / TZ TZ 



lim [C/ Z l tr(* C/^t f/J= lim tr(t E^tf.t C/ Z J7„) = tr It t 
This establishes Eqs. JI2J) of Sect. [3] 



Cndn ■ 
(A6c) 
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B Gaussian target averages 



While direct calculation of specific correlators using the averaging procedure (fTi?|) is in many cases 
straightforward, the calculation can be often simplified significantly by using differential equations. 
This relies on the observation that 

A(... )(y) = .^J^or^^.-.m. (Bl) 

Applied to concrete examples this takes its simplest form if the right-hand side is directly propor- 
tional to (. . .)(Y) itself, but can be useful even in more general cases. We will explicitly address 
the examples needed in Sect. [3j 



Two point projectile 1Z-TZ correlators: Using the notation (|2T|) and a prime to denote a Y 
derivative, straightforward algebra leads to 

d TZ TZ n t ft Ti a TL TZ. t -rO, tl TL H t 

— {tv(U x U y )) = -Q'y^ (tr(i U x U y t ))(Y) = -C n Qy, xy (tv(U x U y )(Y) (B2) 
which is readily solved to obtain 

{tr(U x U y )}(Y) = d n e- c * SY -v . (B3) 
The freedom in the initial condition was used to accommodate the normalization factor dn- 

Three point projectile adjoint-TZ-TZ correlators: These involve several distinct color struc- 
tures. 

:{[U z ] a \v(t U x t U y ))(Y) = -g Y , xz ([fU z } ab tr(i * U x t U y )){Y) 



dY' 1 ZJ v x y> 

+ GY Z y(\i l U z ] ab tv(i Uj U y t ))(Y)+G Y ^ y {[U z ] ao tr(i t U x i U y i ))(Y) 



a l, n H a Tl -rP K t Tli b -R TC a TC* 71 -rP TZ t -R.i 



N c „, x „, / N, 



[On- 



where we have used 

TZ^TZ a 1Z^ TZ a TZ^'JZ^ TZ^ lZ a 7?.* 

t t t = f it + [t , t It 



6 K t 



([t/,] a °tr(t £7 x t (B4) 



TZ^TZ^ lZ a 1Z3 lZ a { J\[ , 

t Cn + ifiakt t = t Cn + ifiak-ifujt =t I C n - -± ) . (B5) 



The nontrivial point is that this holds for any representation TZ. Integrating (|B4p one finds 



.q&TC -R b 7tt 



tr(i 17* t t/„))(y) = c , redRe-^ (ey -" +ey --»- ey -"» ) - CTCSv -»» , (B6) 
again with the free initial condition used to set the normalization properly. 
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C Odderon contributions lead to Casimir scaling violations 



Here we explore the results of applying (|3"T|) to the calculation of dipole and 3-point correlators, as 
was done earlier with the Gaussian truncation. 

It turns out that starting from this level of three-point f-channel correlators generic expressions for 
arbitrary representations 1Z can not be given. This is a consequence of the arbitrarily complicated 
decomposition of a general s-channel TZlZ-dipole into irreducible representations. These start to 
mix in nontrivial ^-dependent ways beyond the Gaussian truncation. 

We have therefore restricted ourselves to 7Z being either the fundamental or the adjoint representa- 
tion. Here it turns out that in (j3"T|) does not contribute at all to correlators of Wilson lines, and 
that G d , as expected, can be thought of as an odderon contribution. This will allow us to compare 
the results we are about to obtain to [54] by counting G d as 0(a s ) (using the quasi-classical count- 
ing) and correspondingly expanding the equations we get to the lowest order in G d , as was done in 
(|38|) . Our results can also be compared to [55] if we keep the G d contributions to all orders. 

For the correlators in question, this contribution generates imaginary parts wherever the representa- 
tion 7Z is not intrinsically real, such as the adjoint representation. For the limited set of correlators 
we are looking at, only the x <-» y antisymmetric combination dY 1 (G Y , yyx — G Y , yxx ) enters. 
For compactness, we will also absorb a constant in the shorthand to be used below. We dchnc 



r 

&Y,xy : ~ ~ J (Gy, yyx — Gy, yxx ) 



(CI) 



with 



N 2 - 4 

C d := -JL— (02) 

characterizing the symmetric "octet" (in SU(3) parlance) in the decomposition of a g_g-dipole into 
invariant multiplcts (see. [69], Sec 9.12 for a systematic treatment that is much more practical than 
most). Note that Q$ xx = 0. 

Let us begin with the three point qqg function. With 1Z the fundamental representation one finds 
= N c Cfe~ [^(S^+S^HS^cv--^)] < 1 — - 



Q° 

^- + n c (g xz 




The expanded expression in the second line serves to recall that higher orders in this expansion in 
powers of Q are beyond the control of our approximation. 

The coincidence limits (|12p , which are at the heart of factorization violations provide relationships 
with dipole correlators also here. Eq. (|C3|) in the limit x = y provides the expression for the 
adjoint correlator 

(tv(U z Ul)}{Y) = ([U z ] ab 2tr{t a U x t b Ul)){Y)= 2N c C f e^ G — (C4) 
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which turns out to be unmodified from the Gaussian truncation and remains completely independent 
of the odderon term Q even without expanding in the odderon contribution. Note that this is 
more stringent than the group theoretical requirement that the gg-dipole has to be real, which 
would have allowed even powers of Q° to appear in an all orders expression in terms of Q° . 

The fundamental correlator, on the other hand, does acquire modifications both to real and imagi- 
nary parts. In accordance with the limits z = x and z = y of (|C3|) one finds 

<tr(l^t))(y) =± { [u x ] ab tr(t a U x t b Ul))(Y) = ±(([U v ] a \r^U x t b U v )) f )(Y) 
W W 

= N c e - c ^y +iS °^y) = N c e~ c ' Q ^ [l-iC, Q? tXy + o (<?° 2 )] . (C5) 

Contrary to what happens in the adjoint representation, this result would emerge from our earlier 
expression for the gg-dipolc using the substitution 

QY, X y -> Gv, X y + i Qy.xy ■ (C6) 



Comparing Eqs. (|C5[) and (|(J4|) we can see that the odderon contribution introduces violation of 
the Casimir scaling of (|23p . We therefore can conjecture that Casimir scaling violating effects are 
due to multiple reggeon exchanges. 

The evolution equation for the qq dipole (Eq. (|lip with 1Z the fundamental representation, after 
insertion of (|C5p and (|C3[) ) leads to 



^(gy^+iG?,^) = J jSz K xxy (l - e ^[(^ + .C) + ( ff , + . 6 f»)-( S .v + .C v )](v)) 



which repeats the structure of (|22|) . again with the simple substitution ()CGI) . 

As we have control only over the terms linear in Q° we expand (|C7[) and use ([22]) to obtain 



(C7) 



^ -^J d2z K **y e-*w-+"--«-.]cn + Gl zy g?. xv ] . (C8) 

Defining the real part of the S'-matrix and the odderon exchange amplitude O by 

S Y , xy = e-°f s *y , Oy. xv = -i C f G? )Xy e- Ci (C9) 

and using these definitions in (|C8|) in the large- A^ c limit (which is needed here just like it was needed 
to derive BK equation (|25p from the GT truncation (|2"2"|) in Sect. 14. ip one derives the non-linear 
evolution equation for the odderon found in [54] (see also [55]) 

-^y OY.xy = ^2 J d 2 zlC XZ y [Oy,xz Sy,zy + Sy,xz Oy.zy — Oy.xy] ■ (CIO) 

In [54,55] the authors discuss how this equation maps onto the BJKP hierarchy [45-48], i.e., onto 
the systematic inclusion of multi-reggeon exchanges in the t-channel. Beyond the low density limit 
[where g° , g (and all higher n-point t-channcl insertions) are small] our procedure provides a 
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generalization consistent with JIMWLK evolution. In the linear regime our solution for O stays 
constant with energy in agreement with [53,54]. 

One of the main properties of the resulting Eqs. (|C8[) and (|C10|) is that Q° = (or, equivalently, 
O = 0) is a stable solution of the equation, and that (as already discussed in [54,55]) non- vanishing 
odderon contributions in the initial condition are erased very quickly due to nonlinear effects. This 
may provide a glimpse of how the Casimir scaling-violating multi-reggeon contributions may be 
erased by non- linear evolution even if they are present in the initial conditions. 

Note that already the inclusion of Q° breaks the Casimir scaling relation (|23|) between the dipolc 
correlators in the fundamental and adjoint representation (see Eqs. (|C5|) and (|C4p V It is only 
natural that higher t-channel n-point exchanges will also contribute to this breaking of Casimir 
scaling. The odderon contribution in our simulations, however, vanishes from the outset: the initial 
conditions necessary to accommodate the total cross section in DIS (as was the case for all our 
simulations) require Q° = 0. 

Still, this discussion does clarify the nature of what we expect to arise as one includes corrections 
to the Gaussian truncation of JIMWLK evolution. 
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